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Abstract 



We prove the existence of a class of perfect-fluid cosmologies which 
have polarised Gowdy symmetry and a Kasner like singularity. These 
solutions of the Einstein equations depend on four free functions of one 
■ space coordinate and are constructed by solving a system of Fuchsian 

00 ' equations. 

o 
in 

Q 1 Introduction 

O 

— . In this paper we use the Fuchsian algorithm to construct a family of perfect- 

q-i fluid cosmologies with polarised Gowdy symmetry and with Kasner like 

^ ■ asymptotics at early times. This family depends on the maximum number of 

bi). free functions for spacetimes within the symmetry class. The technique used 

is to perturb exact Bianchi I solutions in one space-direction and to solve a 
Fuchsian system of equations for the perturbation. The results obtained are 
a generalisation of those in Q in the sense that the symmetry requirement 
has been relaxed a little, but here we require the free data for the field 
equations to be analytic rather than merely C°°. 

The main result is stated as Theorem 6.1 at the end of the paper. 

2 Exact solutions 

The general Bianchi I solution of the Einstein-perfect fluid equations, as 
given in ||] is 

dg 2 = _ B 2(y-l) dT 2 + T 2vr B 2 qidx 2 + ^^^2 + ^^,^2 (1) 

where 

B 2 -~ ( = a + m V~ 7 a > 0, m > (2) 
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Pi + P2 + P3 = 1 , p\ + p\ + P% = 1 , Qi = g - Pi (3) 

(The case a = is an FRW model and is excluded in what follows.) 
The density of the fluid is given by 

4m 2 

When we come to write down Einstein's equations for a metric with po- 
larised Gowdy symmetry we will use conformal coordinates, i.e coordinates 
for which the metric takes the form 

ds 2 = e 2A(t > x 1{-dt 2 + dx 2 ) + R{t, x)e w ^dy 2 + R(t, x)e~ w ^dz 2 (5) 

The model solution (1) may be written in conformal coordinates by making 
the transformation 

t — > t = ( T B pi+ ~f-l(s)s- pi ds (6) 

Then the metric (1) takes the form (5) with 

e A = T Pl B^- pi , e w = T v^B P3 - p2 , R = t^B^ 1 ^ (7) 
and r given implicitly as function of t by the relation (6). 

3 The Einstein-perfect fluid equations 

We will assume that spacetime is filled with a polytropic perfect fluid, so 
that the stress tensor takes the form 

T af3 = {p + P)u a u p + Pg af3 (8) 

with u a u a = -1 and P = (7 - l)p, 1 < 7 < 2. 

With this stress tensor the Einstein evolution equations for the metric 
(5) are 

Ra - R xx = Re 2A p{2 - 7) (9) 
W tt -W xx + R- 1 (R t W t -R x W x ) = (10) 

Att - A xx - \r~\R 2 - Rl) + \(W 2 - W 2 ) = -l 7 e 2A p (11) 
while the constraints read 

R^Rzz - \R- 2 {R 2 t + Rl) - R^R t A t - R~ l R x A x + \{W 2 + W 2 ) 
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= -e 2A p(l+ 7 (v 1 ) 2 ) (12) 

-R-iRtx + R-^RtAx+AtRx^R-tRtRx-lWtWx = -e 2A 1( nPv x (13) 

where v l = e A u l and hence (v ) 2 — (v 1 ) 2 = 1. 

The Euler equations V a T a p = take the following explicit form 

~fv (v°pt + v 1 p x ) + (1 - 7)pt 

+7p{u°(2i; t + A^ 1 ) + ^ - t; 1 ^ + A x v° + ^i; 1 )} 
+ jpv°{v°A t + vl + v 1 A x + R-\v°R t + v 1 R x )} = (14) 

and 

7 vVPt + «Vx) + (1 " 7K + 7P{w°(w* + w°4r)} 
+ 7p{u 1 («? + 2t£ + v°(2A t + iT 1 ^) + ^(A, + iT 1 ^)} = (15) 

4 Inhomogeneous perturbations 

We now seek solutions (p, u 1 , A, it!, W) of (9)-(15) which depend on the max- 
imum number of free analytic functions, namely four, and which approach 
the model forms (l)-(4) for each x as t — > 0. 
To be specific we make the following ansatz 

A = Pl \ogr+ Q -p^j \ogB + tA (16) 

i? = r 1 - pi S^ 1+ l)(l + ^) (17) 
W = c{x) + (p2 - p 3 ) log r + (p 3 - p 2 ) log -B + tT^ (18) 
log p = log n + t (7-i-Pi)/(i-Pi)+^ ( 19 ) 

v 1 = T~ / - 1 - pi (G(x)+t £ 4>) (20) 



where 



£ 2 ~7 = a(x) +m 2 (x)T 2 -~< (21) 



G (x) = ^_ a (p 1 +^)/(2-7) x (22) 

4m 2 7 

1, x , a 3! (37(l-pi) + 2(3pi-l)) f 2(l 
- j (P2 " P3)c, + + (Pi)x (^^yY l0ga + 1 

(23) 

and r(i, x) is determined by the relation (6). 
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a(x) and m{x) are chosen as strictly positive analytic functions, c(x) is 
an analytic function and the Kasner exponent pi(x) is subject to the restric- 
tion pi(x) < 7 — 1 — k for some arbitrarily small constant k. The analytic 
function e(x) is chosen to satisfy 

< e(x) < min ( 7 - 1 - M (24) 
\l-pi 1-pi ) 

We are looking for analytic solutions of the Einstein-perfect fluid equa- 
tions for which A, R, W, 4> and -0 tend to zero as t tends to zero. The 
choice of G{x) ensures that the momentum constraint is satisfied at t = 0. 

It is convenient to write the field equations for A, R, W, 4> and in 
first order form. To do this we introduce the following new variables: 

U = R t , Q = (tA) t , X = R X ,Y = tA x , V = (tW) t , Z = tW x , S = t^R 

In terms of these variables the evolution equations (9)-(ll) take the form 

tX t = tU x (25) 

tY t = tQ x (26) 

tR t = tU (27) 

tA t + A-Q = (28) 

tZ t = tV x (29) 

tW t + W - V = (30) 

tS t + S - U = (31) 

tU t +2U = tX x +2U (l - ^(1 - Pl )B 2 ^ + m 2 + ^\ r 2 ~^ B^^tT^' 1 

+tT Pl - 1 B- ( - Pl+ ^{2X(T 1 - pi B pl+ ^) x + (1 + £)(T 1 - pi Bi ,1+ 3) a . a .} 

+ f^!(2 - 7 )(l + tS)tT 2pi -~<Bl~ 2p ^(eM2tA + ^+(7-i-pi)(i-pi)- 1 ^ ) ) _ !) 
3 (32) 

~\u+ 1 ^ 2 ~ P3 V = ^ + *{(pi logr), + (logSt"^} 
+ i(l+ J R)- 2 tf/ 2 -i[/ (\ - tr^B-^+^il + t^)- 1 ^(1 - Pl )B 2 ^ + m 2 ^pi 
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2m? 



+it r 2(Pi-i)_ B -2(pi+|) (V-pi^pi+Ix + (1 + fl)(r 1_I,1 B Pl+ 5; 
-\tV 2 + hvfa-pz) (^(1 - pi)' 1 - tT^Bl+P 1 -^ + m 2 tT 1+pi -^B-^-P^ 
+\t (c x + ((p 2 - pa) log t) x + ((p 3 - pa) log B)* + 

7B |-2 P1 -7 tT 2 P1 - 7 ^ exp (2ti + ^(T-I-PIKI-Pl)- 1 ^ _ ^ (33) 

tVt + V + ^~ P2, ] u = tZ x + t{c xx + {{P2-Pz) log t) xx + ((P3-P2) log £) ra } 
(1 - Pi) 

f T pi—l R — Pi - 3 i i 

-(c,+(fe-P3)logr) :c +((p3-P2)logi?) :E +Z)(r 1 -^ J B^+3X+(l+ J R)(r 1 -^^+3) ;E ) 



(1 + J2) 

+ F(l - ^-^-^^-^(l - pi)^ 2 - 7 + m 2 r 2 - 7 (pi + 1/3))) (34) 
The Euler equations , after some rearrangement, may be written 

l + 2r 2 V 2 7 t/V 2 V 



1 + 7 r 2 / 3 V 2 



+t 1 - £ r-^ + 7 /^ 2 ^ 2 2 ) ((m 2 ) a; - 7 (T a; + J B a; )+t £+ ' 3 ( 1 -^)"(^+0(/3(l- R )- 1 ) :!; log i )) 

+* 1_e ^ ( 7( y7 7 ^y 2) " 2 ) ((2/3 " ^ iH " Pl ™ 2r " /3 + 0) 
-t^V-" (l + 2r 2 V - ^^ (3«V^ - r 2 V)) ((Pi logrJ^+OogBf-w^+y) 

( ^^2^2 " v °) ((! + + ™V + V3) J B- pi ^r' 3 ) 

+t 1 - £ r^- 2 ^ 2 ^ ^^L - J B- pi -i(l+^)- 1 (r 1 -w^ 1 +^+(l+i?)(r 1 -w J BP 1 +i) :c ) 



(35) 



i7 l + 2r 2 / 3 V> 2 7 

+ 7 t 1+ ^- t T «- 1 ^ f 1 - ^ (Bf-W-7 + m 2 r 2-7 B -|-pi) 

\ 1 + 2r^ V / 
+ tylogt(e + 0(1 - pi)" 1 ). + t 1+ ^" e ((m 2 ) x - 7 (r x + B x )))v ^ 



2 ^ - 1 -2/3,/,2 



+t 1+ ^" (t^{t^ - v°)(( Pl logr) x + (logBf-") x + y) + ^° T ^ " (W0)2) 



1 + i? 



5 Existence and uniqueness of solutions 



(36) 



A careful inspection of the field equations (25)-(36) shows that they may be 
written in the form 

td t u + N(x)u = t 5 H(t,x,u,u x ) (37) 
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where u stands for (U, A,Q, X,Y, R,V, Z, S, (fr,ip), 5 is a strictly positive 
constant and H is continuous in t and analytic in its other arguments. The 
matrix N(x) has positive eigenvalues. 

It follows by d that (37) has a unique analytic solution u with u(0) = 



6 The constraints 

Define constraint quantities Ci, Co by 
Co = R^Rxx - ^R-\R 2 t + R-l) - R- l RtA t - R^R X A X + ~(W t 2 + W 2 X ) 

+ e 2 ^(l+7(^) 2 ) (38) 

Ci = -RT l R tx + R- l {R t A x + A t R x ) + ^R~ 2 R t R x - ^W t W x + e 2A 7 /wV 

(39) 

If the evolution equations (25)-(36) are satisfied then a calculation shows 
that the following hold 

d t c = -a x d - ^c - (40) 

dtd = -d x C - - ^C (41) 

One also calculates that the quantities C$ = tCo, C\ = tC\ tend to zero 
as t tends to zero. These quantities satisfy the following 

td t Co +(jf- 1 )co = -td x C x - (42) 

tdA + (^±-A Cl = -td x C - ^C (43) 

Now R x /R is O(l) and tR t /R = 1 + 0(t s ) for some 8 > 0. It thus follows 
that Co and C\ are identically zero and thus the constraints are satisfied. 
Summarising the results of sections (4)- (6) we have proved the following 
Theorem 6.1 Given two strictly positive analytic functions a(x), m(x), 
an analytic function c{x) and an analytic function pi(x) satisfying 

-g <Pi(x) <7-l-fc (44) 
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for some small k > , there exists a unique solution (g^, p,u a ) of the 
Einstein equations coupled to a 7-law perfect fluid on M 3 x (0, T) satisfying 



ds 2 = e 2A {t,x){-dt 2 + dx 2 ) + R{t,x)e wit ' x) (dy) 2 + R{t,x)e- W{t ' x) (dz) 2 

A = pi log t + log B^ Pl + 0(t) 
R = t 1 -? 1 (1 + O(t)) 

TV = C (x) + (p 2 - P3) log t + (P3 - P2) log B + 0{t) 
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log p = log 



v = 




where r is implicitly given by 
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